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ABSTRACT 

We investigate the physical processes involved in the development of the red sequence 
(RS) of cluster galaxies by using a combination of cosmological -/V-body simulations 
of clusters of galaxies and a semi-analytic model of galaxy formation. Results show 
good agreement between the general trend of the simulated RS and the observed 
colour-magnitude relation (CMR) of early-type galaxies in different magnitude planes. 
However, in many clusters, the most luminous galaxies (Mr ~ My ~ Mx t < —20) 
depart from the linear fit to observed data, as traced by less luminous ones, displaying 
almost constant colours. With the aim of understanding this particular behaviour of 
galaxies in the bright end of the RS, we analyze the dependence with redshift of the 
fraction of stellar mass contributed to each galaxy by different processes, i.e., quiescent 
star formation, and starbursts triggered by disc instability and merger events. We find 
that the evolution of galaxies in the bright end since z ss 2 is mainly driven by 
minor and major dry mergers, while minor and major wet mergers are relevant in 
determining the properties of less luminous galaxies. Since the most luminous galaxies 
have a narrow spread in ages (1.0 x 10 10 yr < t < 1.2 x 10 10 yr), their metallicities are 
the main factor that affects their colours. Their mean iron abundances are close to the 
solar value and have already been reached at z sw 1. This fact is consistent with several 
observational evidences that favour a scenario in which both the slope and scatter of 
the CMR are in place since z ~ 1.2. Galaxies in the bright end reach an upper limit in 
metallicity as a result of the competition of the mass of stars and metals provided by 
the star formation occuring in the galaxies themselves and by the accretion of merging 
satellites. Star formation activity in massive galaxies (stellar mass M* > 1O 1O M0) that 
takes place at low redshifts contribute with stellar components of high metallicity, but 
the fraction of stellar mass contributed since z ft: 1 is negligible with respect to the 
total mass of the galaxy at z = 0. On the other hand, mergers contribute with a 
larger fraction of stellar mass (sa 10 — 20 per cent), but the metallicity of the accreted 
satellites is lower by ~ 0.2 dex than the mean metallicity of galaxies they merge 
with. The effect of dry mergers is to increase the mass of galaxies in the bright end, 
without significantly altering their metallicities. Hence, very luminous galaxies present 
similar colours that are bluer than those expected if recent star formation activity were 
higher, thus giving rise to a break in the RS. These results are found for simulated 
clusters with different virial masses (10 14 — 10 15 /i _1 M Q ), supporting the idea of the 
universality of the CMR in agreement with observational results. 
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It is now well established that galaxies follow a bimodal 
distribution in the colour-magnitude plane, separated into 
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a tight 'red sequence' (RS) and a 'blue cloud'. The RS 
mostly consists of gas-poor galaxies with low levels of star 
formation (SF), prototypically early- type galaxies (ETGs), 
whereas late-type galaxies are typical objects of the blue 
cloud. 

The colour-magnitude relation (CMR) of ETGs can be 
understood as a mass-metallicity relation. The more lu- 
minous (and therefore more massive) galaxies in this re- 
lation have deep potential wells, capable of retaining the 
metals released by supernovae events and stellar winds. 
The CMR seems to be quite universal, since it is followed 
by galaxies in the field as well as in groups and clusters 



(e.g. Visvanathan fc Sandagelll 977; San dage fc Vi svanathan 



19781 ; IReda et al J 12004 lLopez-Cruz et all 12004 iReda et al.. 
20051 : Ide Riicke et alj|2009f ) , but with a larger fraction of red 



galaxies being in denser environments. 

There is controversy on which are the mechanisms re- 
sponsible for the evolution of galaxies to the RS. The phys- 
ical processes invoked can be classified in internal and ex- 
ternal. The former include passive stellar evolution and disc 
instabilities, while the latter involve galaxy-galaxy int erac- 
tions and mergers ([Lin et alfeOlOnRobaina et al.ll2010f ). and 
environmental effects such as fast gravit ational interactions 
('galaxy harassment '. iMoore et al.lll999l ). remo val of the hot 
gas reservoir ('strangulat i on' or 'starvation'. [ Larson et al.l 



1198(1 iBalogh et all I2OO0I ; iKawata fc Mulchaevl |200S)) and 



ram pressure stripping of galacti c gas i Lanzoni et alj|2005t 
iRoedigerj 120081 ; iTecce et alj|2010h . 

Studying the evolution of the slope and scatter of the 
CMR provides constraints o n the SF history of ETGs. 
iBower. Kodama fc Terlevichl (|1998T ) find that the observed 
CMR scatter in local clusters can be reproduced if galax- 
ies are assumed to form the bulk of their stellar content 
at early epochs (z > 1), with a subsequent moderate mass 
growth driven by merger events and residual SF activity. 
This scenario is consistent with hierarchical models of galaxy 
formation. I n this context, using semi-analytical modelling, 
iMenci et al.l (|2008T ) obtain a narrow RS for cluster galax- 
ies, already d efined by z ~ 1.2. At a similar redshift, 
iKavirai et alj (2005) find that the slope of the CMR of clus- 
ter galaxies changes appreciably with respect to its present 
value, although the evolution of the slope in the redshift 
range < 2 < 0.8 is negligible taking into account the er- 
rors. 

Detailed observ ations of the CMR of cluster galaxies 
(see iMei et al.l |2009j, and references therein) show no sig- 
nificant evolution of the zero point, slope and scatter of 
the CMR out to z fa 1.3. In contrast with these results, 
IStott et al.l (|2009r ) find from optical and near-infrared ob- 
servations of cluster galaxies an evolution of the CMR slope 
between z as 0.5 and the present epoch, which they attribute 
to the infall of galaxies into the cluster core that are being 
transformed into RS galaxies. 

Hydrodynamical simulations have also be en used to 
analys e the properties of the CMR at z — 0. ISaro et al.l 
(2006) evaluate the impact of the stellar initial mass func- 
tion (IMF), finding that the Salpeter IMF allows to recover 
both the slope and the normalization of the CMR for galax- 
ies in cluster-sized haloes. Fr om the analys i s of si mulations 
of galaxy groups and clusters, iRomeo et al.1 (2008) conclude 
that the shape of the RS is mainly determined by the spe- 
cific SF in all environments. These authors find that evolv- 



ing galaxies move towards a 'dead' sequence soon after they 
have almost stopped the bulk of their SF. Fainter galaxies in 
clusters keep a significant amount of SF out to very recent 
epochs, and are thus distributed more broadly around the 
RS. 

A linear relation has generally been used to fit the corre- 
lation between luminosity and colour of the CMR, from giant 
to dwarf ellipticals. However, the scatter of these relations 
at lower luminosities, as becomes evident from 



increases 



1997 



observati ons of several cluster s (Coma, Seeker et al.l 
Perseus, IConselice e t al. 200 3]; Forna x, iHilker et al.1 
Karick et alj|2003t iMieske et al.ll2007l ; Hydra, iMisgeld et a] 



2003 



200S| ; Virgo. disker et al.ll200Sl ; Antlia. ISmith Castelli et al 
20081 ). 

Some of these relations seem to be consistent with 
a change of slope from the bright to the faint end. As 
an example, this kind o f fit has already been suggested 
by Ide Vaucouleursl (119611) for t he V irgo c luster, and con- 
firmed later by iFerrarese et al.l (e.g. |2006J). More recently, 
based on Sloan Digi tal Sky Survey (SDSS) imaging data, 
Ijanz fc Lisker] (|2009j) found, also for Virgo, a non-linear re- 
lation that can be described by a 'S' shape over the whole 
range of magnitudes, with the brightest galaxies ( — 21 < 
Mb < —19) havin g an almost constant colour. In the case of 
the Hydra cluster, IMisgeld et al.l (J2008|) fit a linear relation 
to its CMR, but a change of slope is evident for the brightest 
galaxies in that cluster. 

Additional evidence of a tilt towards bluer colours at 
the bright end of the RS arises fr om studies of larg e sam- 
ples of galaxies in the SDSS (e.g. IBaldrv et alj 12004 .120061 ; 
ISkelton et alll2009l '). The results of lSkelton et all (|2009h were 
obtained for the lo cal RS averaged over all environments. 
IBaldrv et al.l (J2006J) examined the dependence of the fit (us- 
ing a 'S '-shaped combination of straight line plus a hyper- 
bolic tangent function) to the mean positions of both the 
RS and blue cloud, and found that unlike the fractions of 
red galaxies, the fits to the RS do not vary strongly with 
environment. 

The detachment of the bright end from the general 
trend denoted by the linear fit to the CMR of cluster galax- 
ies motivates our study. The CMR constitutes one of the 
major tools to test galaxy formation models since galactic 
evolution is evidenced by it. Our aim is to contribute to the 
understanding of the physical processes behind the develop- 
ment of the RS in galaxy clusters, and the special behaviour 
of its bright end. 

Recent observations support strong evolution with cos- 
mic time in the structural properties of the most massive 
(stellar mass M* > 10 Mq) spheroid-like galaxies. Size 
measurements of ellipticals at high redshift (z > 1.5) re- 
veal that these objects are much smaller (by factors ~ 2 
to ~ 4) t han their local co unterparts of comparable stel- 
lar mass (jDaddi et alj 120050 . This result was subsequently 
confirmed by several studies, although the most massive 
ETGs (M* > 2.5 x 10 n M©) seem to have reached their 
internal and dynamical structure earlier and faster than 
lower mass ETGs a t the same redshift, as discussed by 
iMancini et al.l (J2010T >. On the other hand, a mild evolu- 
tion in velocity dispersi on is found for massive galaxies 
IjCenarro fc Tru iillo 2009). On average, the velocity disper- 
sion and surface mass density of massive galaxies at high 
redshift are similar to those of the most dense local ETC 
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l|Cappellari et al.ll2009f ). These observational evidences sug- 
gest that inside-out growth scenarios are plausible, in which 
the compact high redshift galaxie s make up the c entres of 
normal nearby ellipticals (jvan Dokkum et alJl201Cf ). 

From samples of galaxy pairs obtained from different 
surveys, it is found that present-day massive ETGs un- 
dergo, on average , ~ 0.5 major merge rs since z ~ 0.6 — 0.7 
(|Bell et al.ll2006bl ; iRobaina et alj|2010h . These mergers gen- 
erally occur between galaxies which are already on the 
RS, thus involving ETGs containing little cold gas an d 
dust IjWhitaker fc van Dokkuml2008l : Mcintosh et ail2008h . 
These dissipationless (gas-poor, thus called 'dry') major 
mergers are an important channel for the form ation and 
evolu tion of the brightest cluster galaxies (BCGs) l|Liu et al.l 
120091 ). However, the rate of dry major mergers is too low to 
fully explain the formation of spheroidal (and RS) galaxies 
since z ~ 1 (Bund v et al.l 120091) . Minor dry mergers would 
favour the growth of BCGs as inferred from the evolution 
of th eir velocity dispersion-stellar mass relation (|Bernardil 
2009). The minor merger hypothesis also appears to be 
plausible for early-type gala xies at intermediate red shifts 
(0.1 < 2 < 0.8) as shown by iNierenberg et all (|201lh who 
find that ETGs are typically surrounded by several satellite 
galaxies. Minor mergers are also probably responsible fo r the 
size e volution of ma s sive E TGs (|Bezanson et al.l 2009). Be- 
sides, iKavirai et al.l (|201lh demonstrate that the low-level 
of star formation activity in the ETG population at inter- 
mediate redshift is likely to be driven by minor mergers. 

The influence of mergers on the evolution of ETGs in- 
ferred from observational data is in agreement with predic- 
tions obtained from models of g alaxy formation. Usin g semi - 
analytic modelling techniques, iKhochfar fc BurkertJ (|2003l ) 
find that the last major merger of present-day bright el- 
liptical galaxies (Mb < — 21) occurs preferentially between 
bulge-dominated galaxies. Major mergers are responsible for 
doubling th e stellar mass at the ma ssive end of the RS 
since z ~ 1 (Eliche-Moral et al. 2010). On the other hand, 
iBournaud. Jog fc Combea l|2007l ) study the effect of multi- 
ple minor mergers on ETG evolution by using a iV-body 
simulation finding that repeated minor mergers can form 
elliptical galaxies without major mergers events, which are 
less frequent than minor mergers at moderate redshift. Rem- 
nants of these multiple mergers develop global properties 
that depend more on the total mass accreted during mergers 
than on the mass ratio of each merger. This result is also sup- 
porte d by hydrodynamic simulations (e.g. iNaab et alj|2007l . 
I2009T ) which suggest that the growth of massive galaxies 
(M* > 10 1 Mq) by a large series of minor mergers becomes 
more important than the growth via major mergers; minor 
mergers may thus be the main driver for the late evolution 
of sizes and densities of ETGs. Another possible physical 
mechanism responsible for the size evolution of ETGs is the 
expulsion of a substantial fraction of the initial baryo ns, still 
in gaseous form, by quasar activity ( Fan et al.ll2010l) . 

From semi-analytic modelling, IKhochfar fc Silkl (|2006P l 
find that massive galaxies at high redshifts form in gas- 
rich mergers, whereas galaxies of the same mass at low red- 
shifts form from gas-poor mergers. These considerations al- 
low them to reproduce the observed size-redshift evolution 
of massive spheroidal objects. Dissipation associated with 
large gas fractions is the most i mportant facto r determin- 
ing spheroid sizes, as found by iHopkins et al.l (J2009J) us- 



ing high resolution hydrodynamic simul ations. Complemen- 
tarily, also using semi-analytic models, |Pe Lucia fc Blaizoj 
(2007) find that later (z < 0.5) minor dry mergers play an 
important role in the mass assembly of BCGs, which acquire 
half of their mass via accretion of smaller galaxies. 

It is clear that different types of mergers have a strong 
impact on the evolution of ETGs, affecting their morpholog- 
ical and kinematical properties, as well as their star forma- 
tion history. As has been described, many works have been 
devoted to study these aspects but only a few have analysed 
in de tail the role of mergers in the developmen t of the CMR 
(e.g. ISkelton et ail 120091 ; iBernardi et"aHl2010l ). In particu- 
lar, by using a toy model combined wit h galaxy merger trees 
extracted from a semi-analytic model, ISkelton et al.l (|2009f ) 
find that the tilt towards bluer colours of the bright end 
of the RS of local early-type galaxies can be explained by 
the higher fraction of dry mergers suffered by bright galax- 
ies with respect to fainter ones. Galaxy colour is not ex- 
pected to change during dry mergers, since there is no asso- 
ciated SF. Galaxies then move along the CMR as the mass 
of the system increa ses, with their colour remaining fixed 
l|Bernardi et alj|2007f ). 

In this work, we investigate this issue in detail us- 
ing the semi-analytic model of galaxy formation and 
evolution SAG (acronym for Semi- Analytic Galaxies; 
lLagos. Cora fc Padillal |200S| . hereafter LCP08), which is 
combined with hy drodynamica l cosm ological simulations of 
galaxy clusters by iDolag et al.l (|2005T >. We focus our study 
on galaxies lying on the RS, which are selected by a colour 
criterion. The samples thus obtained mostly include early 
type galaxies, making our conclusions valid for the under- 
standing of the development of the CMR of ETGs. Thus, we 
hereafter use the term RS to refer to colour-selected galaxies 
from our simulations, and the term CMR to refer to observed 
relations for ETGs. 

The tilt in the bright end of the RS emerges naturally 
from the model once it has been calibrated to recover nu- 
merous observed galaxy properties, both locally and at high 
redshift. We quantify the contribution to the mass and metal 
content of galaxies located in the RS given by quiescent SF 
from the cold gas available in the galaxy disc, starbursts trig- 
gered by disc instability and merger events, and the stellar 
mass accreted from satellite galaxies during mergers, which 
are identified as minor/major and wet/dry ones. 

This paper is organised as follows. In Section [2] we 
briefly describe the semi-analytic model used in this work 
and give details of the cluster simulations. Section[3]presents 
the comparison between the simulated RS and observed 
CMR followed by early-type cluster galaxies in different 
magnitude planes, and the analysis of the luminosity- 
metallicity relation followed by galaxies in the RS. Section [4] 
describes the formalism used to quantify the contribution 
of different processes to the mass and metallicity of galax- 
ies; the corresponding results are presented and discussed. 
Finally, in Section [5] we summarize our main findings. 



2 HYBRID MODEL OF GALAXY 
FORMATION AND EVOLUTION 

In this work, we focus our analysis on the RS of cluster galax- 
ies. We investigate the development of the RS by applying 
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a numerical technique which combines cosmological non- 
radiative Af-Body/SPH (Smoothed Particle Hydrodynam- 
ics) simulations of galaxy clusters in a concordance A Cold 
Dark Matter universe and a semi-analytic model of galaxy 
formation and evolution. In this hybrid model, the outputs 
of the cosmological simulation are used to construct merger 
histories of dark matter (DM) haloes and their embedded 
substructures, which are then used by the semi-analytic code 
to generate the galaxy population. In this Section we briefly 
describe the hybrid model. 



2.1 Simulated galaxy clusters 

We consider two sets of simulated galaxy clusters, C14 and 
C15, which contain clusters with virial masses in the ranges 
~ (l.l-1.2)xlO 14 /i _1 M and~ (1.3-2.3) x 10 15 hT x M , 
respectively. The clusters in the C14 set are named gl542, 
g3344, g6212 , g676 and g914,w hereas those in C15 are gl, 
g8, g51 (see iDolag et al.l 120051 . for details). These clusters 
have been initially selected from a DM-only simu lation with 
a box size of 480 hT 1 Mpc (|Yoshida et all 12001] ) . for a cos- 
mological model characterized by Q m = 0.3, Oa = 0.7, 
H = 70 km s" 1 Mpc" 1 , fi b = 0.039 for the baryon den- 
sity parameter, and a$ = 0.9 for the normalization of the 
power spectrum. The Lagrangian regions surrounding the 
selected clusters have been re-simulated at higher mass res- 
olutio n by applying the 'Zoomed Initial Condition' tech- 
nique (iTormen et al. Il997h and the Tree-SPH GADGET-2 code 
|Springel et al.l 2005). The mass resolution is the same for 
all simulations, with masses of DM and gas particles m<im = 
1.13 x W 9 h- 1 M e and m gas = 1.69 x 1O 8 /i _1 M , respec- 
tively. As for the force resolution, the Plummer-equivalent 
gravitational softening is fixed at s = 5/i _1 kpc in physical 
units at redshift z ^ 5, while it switches to comoving units 
at higher redshifts taking a value of e = 30 h _1 kpc. 

Although the simulations are hydrodynamical, we only 
use the kinematical information provided by DM particles 
in order to identify DM haloes. To this aim, we consider 
that DM particles have their masses increased to its orig- 
inal value, since gas was introduced in the high-resolution 
region by splitting each parent particle into a gas and a 
DM particle. Dark haloes are first identified as virialized 
particle groups by a friends-of-friends (FO F) algorithm. The 
SUBFIND algorithm |Springel et alj|200ll ) is then applied to 
these groups in order to find self-bound DM substructures, 
which we call subhaloes. Subhaloes are the sites of galaxy 
formation, and their merger trees are built by extracting 
from the simulations all subhaloes with 10 or more bound 
DM particles, since sm aller ones are dynamically unstable 
JKauffmann et alJlJ999h . 



2.2 Semi-analytic model SAG 

The semi- analytic model SAG used here is based on that de- 
scribed bv lLagos et al.l ((2008). The physical processes con- 
sidered in this model are the cooling of hot gas as a result 
of radiative losses, quiescent star formation, and starbursts 
during disc instability events and galaxy mergers. Feedback 
from supernova explosions and active galactic nuclei (AGN) 
are also considered. The mass of hot gas available in a halo 
is determined at each simulation output on which the semi- 



analytic model is applied. This mass is defined as the bary- 
onic fraction of the virial mass of the halo minus the total 
mass of baryons in the form of cold gas and stars and the 
mass of black holes already present within the halo. The 
model tracks the circulation of metals between the differ- 
ent baryonic components, that is, hot diffuse gas, cold gas 
and stars, taking int o account t he mass evolution of different 
chemical elements (jCoral|2006l ). The version of SAG used in 
this work incorporates the mo re detailed c a lculat ion of disc 
galaxy scale radii described in lTecce et a.1.1 (| 20101 ) ■ 

After tuning the free parameters involved in the dif- 
ferent physical processes considered, the model is able to 
reproduce several observational constraints simultaneously, 
such as galaxy luminosity functions, relations between cen- 
tral black hole mass and other properties of the host bulge, 
quasar luminosity function, stellar mass functions at differ- 
ent redshifts, and mass-metallicity, colour-ma gnitude and 
disc size-luminosity relations (see LCP08 and iTecce et all 
l20ld for details). 

During their evolution, galaxies can acquire stellar mass 
by SF, either in quiescent form by converting their cold 
gas into stars, or in starbursts which can be triggered by 
disc instability events or galaxy mergers. During mergers, a 
galaxy's stellar mass also increases due to the stars accreted 
from the merging satellite. In SAG, as in many other semi- 
analytic codes, only central galaxies have cooling flows; the 
cold gas accreted is assumed to settle into the galactic disc. 
A galaxy's mass can subsequently be reduced through recy- 
cling due to stellar winds and supernovae explosions. The 
latter also reduce the amount of cold gas available for fur- 
ther SF through energy feedback. AGN activity also injects 
energy into the hot phase and decreases the flow of cooling 
gas into the galaxy, thus regulating SF as well. 

The complex combination of all these processes is re- 
flected in the resulting metallicities and luminosities of 
galaxies. In the present study, we track the increase of galaxy 
mass considering, separately, the contribution of new stars 
formed during SF episodes of different types (quiescent mode 
and starburst modes) and stars accreted from satellite galax- 
ies during mergers. SAG considers both major and minor 
mergers, and in this work we also identify them as either 
wet or dry mergers. In the following, we give a short de- 
scription of the parameters that characterize different types 
of mergers, and refer the reader to LCP08 for details on 
quiescent SF and disc instability events. 



2.3 Galaxy mergers and starbursts 

In a hierarchical scenario of structure formation, mergers of 
galaxies are a natural consequence of the mergers of DM 
haloes in which they reside, and play an important role 
in determining the mass and morphology of galaxies. The 
galaxy catalogue is built up by applying the semi-analytic 
model to the detailed DM subhalo merger trees extracted 
from the cluster simulations. After the identification of DM 
substructures within FOF groups, the largest subhalo in 
a FOF group is assumed to host the central galaxy of the 
group, located at the position of the most bound particle of 
the subhalo. These galaxies are designated as central galax- 
ies. Central galaxies of other smaller subhaloes contained 
within the same FOF group are referred to as halo galax- 
ies. The subhaloes hosting these galaxies are still intact after 
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falling into larger structures. There is a third group of galax- 
ies generated when two subhaloes merge and the galaxy of 
the smaller one becomes a satellite of the remnant subhalo. 
We assume that these satellite galaxies merge with their cor- 
responding subhalo central galaxy on a dynamical friction 
time-scale. 

Mergers are classified as major or minor according to 
the ratio between the baryonic mass (stellar plus cold gas 
mass) of the accreted satellite galaxy and of the central 



galaxy, /„ 



= M sat /M c 



A major merger occurs 



when /merge > 0.3. Otherwise, we are in the presence of 
a minor merger. 

In the case of a major merger, all stars present in the 
merging galaxies are rearranged into a spheroid. Addition- 
ally, all the cold gas available in the remnant is consumed 
in a starburst as a result of the perturbation introduced by 
the merging satellite. This perturbation drives the cold gas 
into the bulge component of the remnant, where it is com- 
pletely transformed into stars. Thus, major mergers lead to 
the creation of an elliptical galaxy. 

On the other hand, during minor mergers all the stars 
of the merging satellite galaxy are added to the spheroid of 
the remnant, whereas the stellar disc of the accreting galaxy 
remains unaltered. The presence of a starburst during a mi- 
nor merger depends on the gas mass fraction of the disc 



of the central galaxy, / c Xcja 



AI, 



ColdGas 



/MIS 



where 



Mdisc = Mstellar - M B ulgc + M C oldGas , M S tcllar IS the to- 
tal stellar mass of the galaxy, MB„i ge is the mass of the 
bulge (formed only by stars), and McoidGas is the cold gas 
content of the galaxy. If /coHgL > /gas.minor, all the cold 
gas in the merging galaxies is assumed to undergo a star- 
burst; stars formed in this fashion are then added to the 
bulge of the remnant. The inclusion of the minor merger 
starbu rst threshold (,f ga a. m i nor) is motivated by the sugges- 
tion of Irlernauist fc Mihos (1995) that gas-rich discs should 
be suscep tible to bursts trigge red by the accretion of small 
satellites. iMalbon et al.l (|2007f ) consider this SF channel in 
a semi-analytic code; varying the minor merger starburst 
threshold with respect to the value adopted by them, we 
find that we are able to reproduce the observed luminosity 
function with / gaSl minor = 0.6 (as done in LCP08). This type 
of merger is considered a minor wet merger since it leads to 
SF. In this case, the stellar mass of the central galaxy is in- 
creased by the new stars formed and by the stellar content 
of the accreted satellite galaxy. 

On the contrary, if / C omgL < /gas,minor then no burst 
occurs, and we classify it as a minor dry merger. Further- 
more, bursts do not occur if the satellite is much less massive 
than the central galaxy (/merge < /burst, with / bu rst=0.05), 
regardless of the amount of gas in the central disc. Since 
this type of minor merger does not trigger starbursts, it is 
also considered a minor dry merger. The choice of these pa- 
rameters allows us to reconcile observational data with the 
properties of the galactic populations given by SAG. 

For the sake of our study, we also need to classify major 
mergers as either wet or dry. This classification does not af- 
fect how SAG works, and is only of interest for the identifica- 
tion of different types of processes that contribute to galaxy 
formation. Then, once we have identified a major merger, we 
classify it as wet or dry depending on the gas fraction of the 
disc of the remnant. Then, if /cokiGas < /gas, major the major 
merger is dry, and it is wet otherwise. We adopt the value 



/gas, major = 0.4; further discussion regarding the choice of 
this parameter will be given in Section [4] 



3 COLOUR-MAGNITUDE RELATION 

On the basis of the results obtained by applying our semi- 
analytic code SAG to the sets of simulated galaxy clusters 
C14 and C15, we are able to construct the corresponding 
RS, whose bright end is mostly populated by ETGs. In the 
simulations, galaxies belonging to the RS are selected ac- 
cording to the empirical redshi ft-dependent sepa ration be- 
tween colour sequences given bv lBell et al.l l|2004r ). based on 
the bimodal distribution of observed galaxy colours out to 
Z ss 1. In this way, galaxies redder than 

(U-V) = 1.15 - 0.31 z - 0.08 (A/ v - 51ogh + 20), (1) 

for < z < 1, are considered as belonging to the RS. In 
the following, we adopt z = to obtain simulated RS at the 
present epoch and compare them with observed CMRs of 
ETGs. As we will be basically dealing with the CMR slope 
but not with its zero-point, we do not need to convert obser- 
vational data from absolute into apparent magnitudes and 
vice-versa. In the model, neither magnitudes are influenced 
by extinction nor the colours are affected by reddening. 



3.1 Slopes of the colour-magnitude relations in 
different bands 

Our semi-amalytic model provides the magnitudes of galax- 
ies in the Johnson-Morgan system. In order to compare with 
observed data, they have to be transformed to other photo- 
metric systems. Fig. Q] shows the RS obtained for one of the 
more massive simulated clusters (C15_3) in three different 
colour-magnitude planes: R vs. U — R, R vs. (B — R), and Ti 
vs. (C — Ti). We also show the RS of a less massive cluster 
(C14_2) in the Ti vs. (C - Ti) plane in the right bottom 
panel. The R magnitude is in the Cousins system, while the 
C and Ti magnitudes correspond to the Washington photo- 
metric system. The (C — Ti) integrated colour has proved to 
be quite sensitive to metallicity for o ld, single stellar popu- 
lations like globu lar clusters (e.g., see lHarris fc Harris 2002; 
iForte et al.ll2007l ). while the (U — R) integrated colour is age 
sensitive. Besides, the RS in the V vs. (V — I) plane is also 
considered in the analysis, but not shown in this figure. 

The R and / magnitudes in the Johnson-Morgan sys- 
tem are transfor med into the Cousins system through the 
relations given bv lFukugita et all (|l995T ) for ETGs. The pho- 
tometric conversion from Johnson-Morgan V minus Cousins 
I colours into the Washington (C — Ti) ones is performed 
through the transformation (V - I) = 0.49(C - Ti) + 0.32, 
given bv lForbes fc Forte! (J200l|) for globular clusters, under 
the assumption that ETGs are mainly old stellar systems as 
well. Finally, Ti magnitudes are obtained from the Cousins 
R ones by ap plying the relation R — Ti ~ —0.02, given by 
iGeislerl (|l996h . 

It can be clearly seen from Fig. [T]that the simulated RS 
of the C15 cluster, in the three colour-magnitude planes, 
present a break with brighter galaxies detached from the 
general trend denoted by the fainter ones. This break oc- 
curs at approximately the same magnitude in the different 
systems (M]? rcak ~ M^ rcak ~ M^ eak « -20; galaxy mass 
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Figure 1. RS obtained from the semi-analytic model SAG applied to the C15_3 and C14_2 simulated galaxy clusters (see Table 1) in 
different photometric systems: R vs. (U — R) and R vs. (B — R) for C15_3 cluster (top), and Ti vs. (C — Ti) for C15_3 cluster (bottom 
left) and C14_2 cluster (bottom right). Symbols are coloured according to the iron abundance of the stellar component of the simulated 
galaxies. Least-squares fits to the faint end of the simulated RS and to the whole relation are represented by blue solid lines and black 
dashed lines, respectively. The fit to the bright end of the simulated RS of the C15 cluster is also shown in the plane T\ vs. (C — Ti)(red 
dotted line). 



of ~ 10 10 Mq), being more evident in the Washington sys- 
tem, with brighter galaxies displaying an almost constant 
colour (C — Ti ~ 1.5). Taking such break as a reference, we 
will refer hereafter to the galaxies brighter/fainter than the 
magnitude of the break as the 'bright/faint end' of the RS. 

Two fits are performed to the simulated RS in each pho- 
tometric system, considering the whole relation and the faint 
end only. These fits are shown in Fig.[T]for the C15_3 cluster 
and for C14_2 cluster as black dashed and blue solid lines, re- 
spectively. Table [l] presents the slopes of these least-squares 
fits performed to the RS of the eight simulated clusters de- 
scribed in Subsection 12. II Note that the slopes depicted in 
Table[T] designated as b in the following analysis, correspond 
to RS where the colour is the independent variable. In order 
to compare with those obtained from observed CMRs, we 
will have to invert the 'observed' slopes (1/6) in the cases 



that CMRs are calculated taking the colour as the depen- 
dent variable (i.e. as if the diagrams displayed in Fig.[l]were 
rotated 90°). 

With regard to t he observed CMRs, a wide cluster 
sample is presented bv lLopez-Cruz et al.l (12004! ) ■ who study 
the CMRs of 57 low-redshift cluster galaxies, with differ- 
ent richness, cluster types and X-ray luminosities. The fits 
to the CMRs are performed in the (B — R) vs. R plane 
(see their table 1), from which we calculate an average slope 
1/&BR = -0.051 ± 0.002, or 6 B r = -19.61 to compare with 
our simulations. This falls within the range of values pre- 
sented in Table [T] for the corresponding 'all CMR' fits. 

iMisgeld et all l|200ct) fit the CMR of the early-type 
dwarf galaxy sample (Mv > —17) of Hydra I cluster in 
the (V — /) vs. V plane, and get a slope l/6vi = —0.039 
(rms error of the fit equal to 0.12), or 6vi = —25.64. They 
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Table 1. Slopes calculated from least-squares fits to both the faint end of the RS and the whole relation, for the different simulated 
clusters in several photometric bands. 



R vs. B - R 



V vs. V ■ 



Ti vs. C - Ti 



R vs. U - R 



Cluste 



Faint end 



All 



Faint end 



All 



Faint end 



All 



Faint end 



C14.1 (gl542) 
C14.2 (g3344) 
C14.3 (g6212) 
C14.4 (g676) 
C14.5 (g914) 



-16.92 
-18.55 
-19.23 
-17.43 
-16.39 



-19.90 
-20.00 
-20.88 
-20.38 
-18.99 



-20.12 
-19.97 
-24.65 
-26.47 
-21.48 



-25.22 
-27.85 
-30.37 
-28.95 
-25.31 



-10.79 
-12.24 
-12.36 
-11.45 
-10.47 



-12.29 
-13.77 
-13.80 
-13.83 
-12.35 



-8.55 
-12.24 
-9.58 
-8.75 
-8.30 



All 



-9.85 
-13.77 
-10.07 
-9.79 
-9.25 



C15.1 (gl) 


-18.10 


-21.36 


-24.78 


-29.72 


-11.48 


-13.96 


-9.00 


-10.32 


C15.2 (g8) 


-18.11 


-20.69 


-24.12 


-28.84 


-11.55 


-13.71 


-9.15 


-10.21 


C15.3 (g51) 


-18.96 


-21.78 


-25.57 


-31.19 


-12.38 


-14.62 


-9.33 


-10.48 



also estimate the slope of the CMR for Local Group dwarf 
ellipticals (dE) and dwarf spheroidals (dSph) in the same 
colour-magnitude plane as 1/fevi = —0.038 (rms of the fit 
0.09), or bvi = —26.32. Still in the same plane, a slope 
l/b V i = -0 . 033 ± 0.004, or 6yi = -30.30, is given by 
iMieske et ail l|2007l ) for the dE galaxy sa mple (My > -17) 
in the Fornax cluster. JMisgeld et all (|2009l ) analyse the CMR 
followed by dwarf ETGs in the Centaurus cluster and obtain 
1/fevi = -0.042 (rms of the fit 0.10), or bvi = -23.81. Ex- 
cept for the slightly lower value from Mieske et al., all the 
slopes are within the range covered by those of the 'faint 
end' fits listed in Tabled] 

The Perseus cluster is studied bv lde Riicke et al.l (| 20090 
with HST/ACS data; they present a global CMR that in- 
cludes several other groups and clusters, obtaining a com- 
mon fit with slope l/&vi = -0.033 ±0.004, or &vi = -30.30, 
which is also in agreement with 'all CMR' values from Table 

m 



ISmith Castelli et all (|2008l ) obtained the CMR of ETGs 
located in the central region of the Antlia cluster, using the 
Washington photometric system. This CMR is characterised 
by the least-squares fit Ti = (38.4 ± 1.8) - (13.6 ± 1.0) (C - 
Ti). This latter fit corresponds to the whole galaxy range, as 
a break is not distinguishable in this observed CMR. Note 
that this is so because the brightest galaxies in this sample 
hardly reach the magnitudes above the break detected in 
the simulations. 

The linear fit to the bright end of the RS in this colour- 
magnitude plane for the C15 cluster is represented in the 
bottom left panel of Fig. Q] by a red dotted line, and gives 
M Tl = 83.65 - 67.92(C - Ti). Its slope clearly differs from 
the one corresponding to the faint end, which is fit by the 
relation M Tl = 0.82 - 12.38(C - Ti), shown with a blue 
solid line. This latter fit is closer to the one applied to the 
whole RS, for which we obtain M Tl = 2.19 - 14.62(C - Ti). 
The slope of the fit to the whole RS of the C14_2 cluster 
is —13.77 (see Table [TJ in very good agreement with the 
observed value —13.6 of the slope for Antlia. Similar re- 
sults have been found in all the simulated clusters, as can 
be seen from Table [l] This supports the idea of the univer- 
sality of the CMR as in f erred from observational results (e.g 
lLopez-Cruz et alj|2003 : lde Riicke et al.ll2009h . 

The detachment of the bright end (-24 < M Tl < -20) 
with respect to the linear fit to the faint end is clear for the 
three massive simulated clusters (C15), although the change 
of slope in the RS of the less massive clusters (C14) is not so 



evident. This is mainly due to the lack of objects in the high 
luminosity range of C14 clusters and the intrinsic scatter of 
the relation. Only two small groups (C14_l, C14_2) show a 
slightly tighter relation at Mt ± ^ — 20, where the change of 
slope becomes perceptible. As mentioned above, this effect 
is present in some observed CMRs as well. For instance, it 
can be seen in the CMR of Hydra I cluster presented by 
iMisgeld et"afl |200ct ) (their Figs. 2 and 10), though it is not 
particular ly noticed by the aut hors. The studies of the Virgo 
cluster bv ljanz fc Lisk cr ( 2009) and of a large sample of local 
ETGs bv lskelton et all l J200g ). both based on SDSS data, 



clearly show a tilt towards bluer colours at the bright end 
of their CMRs. 

We find a very good agreement between the slopes of 
the observed CMR and the simulated RS, specially for galax- 
ies fainter than the break M^' * ~ —20. It is remarkable 
that this global agreement, as well as the tilt of the bright 
end towards bluer colours, emerge naturally from SAG after 
calibrating the code to satisfy several other observational 
constraints simultaneously, as mentioned in Section [2] 



3.2 Luminosity-metallicity relation 

One of the main properties also affected by the physical 
processes driving the development of the CMR is the galaxy 
metallicity. Simulated galaxies that form the RS shown in 
Fig. [T] are identified by different coloured symbols according 
to the metallicity range in which they lie. The metallicity 
of each galaxy is characterised by the iron abundance of its 
stellar component by using the index [Fe/H], where we are 
considering the solar val ue of iron abundanc e by number 
(Fe/H) Q = 2.82 x 10" 5 (JAsplund et al.ll2005l ). We can see 
that galaxies at the bright end are the most metal rich ones 
with [Fe/H] > —0.25, some of them reaching values as high 
as [Fe/H] ~ 0.43; the mean value of their metallicity distri- 
bution is ([Fe/H]) = —0.03 ± 0.05. As we move down along 
the RS, galaxies become fainter, bluer and less chemically 
enriched. This plot clearly reflects the fact that the RS can 
be interpreted as a mass-metallicity relation. 

We can further test whether a luminosity-metallicity re- 
lation is present. In the upper panel of Fig. [2] we show this 
relation for the same simulated galaxies that form the RS 
of the C15_3 cluster shown in Fig. Q] In order to compare 
our results with observed data, we include [Fe/H] abun- 
dances of early-type galaxies belonging to the NGC 5044 
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Figure 2. Top panel: Luminosity-metallicity relation for the 
same simulated galaxies in the RS of the C15_3 cluster shown 
in Fig. [l] Differently coloured symbols identify galaxies within 
different ranges of age; black open diamonds represent the [Fe/H] 
abundanc es of early-type gal axies belonging to the sample de- 
scribed in Trag er et al] (120001 1 - whereas bla ck open squares cor - 
respond to galaxies in the NGC 5044 group JMendel et al.ll2009r) . 
Bottom panel: RS of simulated galaxies in the magnitude plane 
T\ vs. (C — Ti) colour-coded according to the range of ages in 
which they lie. 



group | Mendel et al] 12003 ) , and of a sample described in 
llrager et al] ( 20001 ). The latter contains cluster members 
(six of the Virgo cluster, eleven of the Fornax clusters, one 
of the rich cluster Abell 194), but most of the galaxies are 
in small groups of varying richness. T he metallicity of thi s 
set of galaxies were re-determined bv iTraeer et al] (2008); 
[Fe/H] abundances are then obtained fr om these metallicity 
values using a transformation given bv lTrager et al] (|2000r ) 
(see their Eq. 1). In all cases, we have transformed the ab- 
solute magnitudes to the T\ band through the relations of 
iFukugita et al] | | 1995ft . As can be seen, the observed values, 
which cover a wide range of magnitudes and abundances, 
overlap the simulated relation. The agreement is particularly 



good for galaxies in the faint end of the RS (Mt x > —20), 
although we note that the observed abundances of some 
galaxies in the bright end are slightly larger than those ob- 
tained from our model. This discrepancy may be due to the 
fact that most of the galaxies in the observed sample are 
in relatively low-density environments, having different SF 
and enrichment histories than those residing in massive clus- 
ters. The fact that the metallicities of galaxies at z — sat- 
isfy observational constraints is particularly relevant for our 
study. It supports the use of the chemical history of galaxies 
as a tool to help understand the development of the RS and 
its special feature at the bright end. 

Galaxies in the luminosity-metallicity relation of Fig. 
[2] have been colour-coded according to the range of ages 
in which they lie. Ages have been estimated considering 
the stellar mass weighted mean of the birth-time of each 
single stellar population. A clear correlation between age 
and metallicity is evident for the least luminous galaxies 
(Mt ± > —16). Younger galaxies are more metal poor, lying 
on the left side of the luminosity-metallicity relation, while 
older galaxies have higher iron abundances being located to- 
wards the right. As we move to brighter galaxies along this 
relation, ages and metallicity be come anticorre l ated, con- 
sistent with the trend found by iGallazzi et al] (|2006T ) for 
galaxies in the SDSS at fixed velocity dispersion. This an- 
ticorrelation might explain the modest scatter of the bright 
end of the CMR (|Trager et al]|2000| ). which is characterized 
by a negligible spr ead in age, al s o con sistent with obser- 
vational results by IGallazzi et al] (|2006t ) . Most galaxies in 
the bright end of our simulated RS share very similar ages 
(1.0 x 10 10 yr < t < 1.2 X 10 10 yr). This can be seen in the 
lower panel of Fig. [2] where galaxies along the RS in the 
Ti vs. (C — Ti) plane are identified with different coloured 
symbols according to their age. From this figure, it is clear 
that, in the faint end, the scatter in colour about the RS 
for a fixed l uminosity is cau s ed ent irely by age variations, as 
assumed bv lBernardi et al.1 ([2005), where younger galaxies 
correlate with bluer colours. As we move to higher lumi- 
nosities, galaxies of a given age become redder because they 
are more metal rich. The effect of age differences on the 
final colours of galaxies in th e bright end of the RS is com- 
pletely negligible (see fig. 1 of lBruzual fc Cha riot 2003]), and 
we can attribute to metallicity effects the behaviour of the 
galaxies in the bright-end of the RS. In the following, we an- 
alyze the physical reasons that make the old galaxies in the 
bright end reach an upper limit in metallicity. Our aim is 
to demonstrate that this metallicity effect is the responsible 
of the colours achieved by massive galaxies which are bluer 
than those that would be obtained from an extrapolation of 
the faint end of the RS. 



4 PHYSICAL PROCESSES INVOLVED IN THE 
DEVELOPMENT OF THE RS 

The galactic properties that determine the RS, that is, 
masses and metallicities, are the result of a complex com- 
bination of different physical processes involved in the for- 
mation and evolution of galaxies. Thus, in order to identify 
which mechanisms are responsible for the different features 
of the RS, and in particular the break at the bright end, 
we track the evolution of the masses and metallicities of 
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Figure 3. Stellar mass fractions of 'stars' and 'sat' contributions 
reached at z = by galaxies in the C15 clusters for the different 
magnitude bins. 
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of the metallicities of galaxies, given by the ratio Fe/H. For 
each process considered, we keep track of the mass of Fe 
and H received by a galaxy from formed and/or accreted 
stars at a given redshift as a result of a certain process; 
thus, we have the quantities Fe^° ccomp and H^° ccomp , re- 
spectively. In order to achieve a better visualization of the 
mechanisms that determine the properties of galaxies along 
the RS, we divide the simulated RS into six bins of one mag- 
nitude from Mt x — —16 to Mt x — —22, and a seventh one 
from Mt x — —22 to Mt x — —24. This last bin comprises 
two magnitudes to improve the statistic analysis, in view of 
the low number of very luminous galaxies in all simulated 
clusters. Thus, bin = 1 and bin = 7 correspond to the mag- 
nitude ranges -24 < M Tl < -22 and -17 < M Tl < -16, 
respectively. We analyze the evolution with redshift of the 
fractions of stellar mass contributed by different processes 
to galaxies in a given magnitude bin, and of their mean 
metallicity. To this aim, we define several quantities that 
are included in these calculations. 

The total mass of stars acquired by galaxies in a mag- 
nitude bin b at redshift z arising from a given component 
contributed by a particular process is: 



Nb.z 



q i j-proc, comp 



= E SM < 



proc, comp 

9,z i 



(2) 



9=1 



the stars added to each galaxy by different processes: quies- 
cent star formation from the cold gas available in the galaxy 
disc, starbursts during mergers and disc instability events, 
and the stellar mass accreted from satellite galaxies during 
mergers. 

In the following, the stellar mass generated during any 
type of star formation event is denoted as 'stars', while the 
stellar mass accreted during mergers, already present in the 
satellite galaxy, is referred to as 'sat'. Thus, 'stars' and 'sat' 
indicate the class of stellar 'component' contributed by a 
given 'process'. Hence, the stellar mass of a galaxy g at a 
certain redshift z that results from the addition of a certain 
stellar component is denoted by SM$™ C ' comp . 

Therefore, we consider a set of variables that take into 
account the different contributions of stellar mass to a given 
galaxy, according to the process involved. Specifically, 

(i) stellar mass from quiescent star formation: 



• SMZ 1 



accent, stars 



(ii) stellar mass from disc instability events: 

• SM^f' stars 
(iii) stellar mass from merger events: 



• SMT™ 1 dry - sat 

n i i-minor wet, sat 
q ifminor wet, stars 






or dry, sat 
or dry, stars 
or wet, sat 
or wet, stars 



With all these separate variables, we track the mass assem- 
bly history of galaxies and the evolution with redshift of 
the stellar mass fractions contributed by the different pro- 
cesses. Similar variables are defined to follow the evolution 



where Nb, z is the number of galaxies within each bin b at 
each snapshot of the simulation characterized by a redshift 
z. Then, considering the combined contribution of the two 
stellar components ('stars' and 'sat') by any given process, 
we have 



SM proc = £ M Proc, stars + £ M PK>c, sat _ 



(•3) 



Finally, the total stellar mass acquired by all galaxies in a 
given magnitude bin that results from the contribution of 
all the involved processes is 



(i)+(ii)+(iii) 



SM, 



where iV D 



E 

proc=l 



SM» r z oc , 



(4) 



is the total number of processes considered 
grouped in the sets (i), (ii) and (iii), as described previously. 
In the following, all results shown correspond to aver- 
ages over the three C15 clusters; error bars are not included 
in the plots because they are extremely small. We have also 
evaluated the situation considering C14 clusters. 



4.1 Evolution of stellar mass fractions 

Taking into account the definitions given above, we estimate 
the fraction of stellar mass accreted by galaxies in each mag- 
nitude bin at a given z resulting from the contribution of 
each of the processes considered, distinguishing between the 
two different components, as 



SMFrac pr z oc ' comp 



q ii j-proc, comp 



S< +C 



i) + (iii)' 



(5) 



These stellar mass fractions can be grouped according to the 
aspect we want to analyse. 

We first consider the fractions of the stellar mass given 
by new formed stars, on one hand, and accreted stars, on the 
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Figure 4. Evolution of the mean stellar mass fractions contributed by different processes for galaxies within magnitude bins 1, 3 
and 6, normalized to the total contribution to the stellar mass at each redshift the three C15 clusters are considered. Left panel: 
galaxies with manitudes within the range —18 < Mt, ^ — 17, magnitude bin 6. Middle panel: galaxies with manitudes within the range 



—21 < Mt x ^ —20, magnitude bin 3 (critical bin). Right panel: ; 
bin 1. 



laxies with magnitudes within the range — 24 < M^ ^ — 22, magnitude 



other, regardless of the processes that contribute to them. 
Thus we have 



SMFrac^° z mp = ^ 



SMFrac pr z oc ' comp , 



(6) 



proc=l 



for the 'stars' and 'sat' components. The dependence of these 
fractions with the magnitude bins, in which galaxies along 
the RS lie, is presented in Fig.[3]for the two components, con- 
sidering the values of stellar mass reached at z = 0. As can 
be seen, the stellar mass received by all galaxies along the RS 
is mainly composed by the 'stars' contribution ranging from 
~ 95% for the least luminous galaxies (Mt x ^ — 17, lying in 
bin 7), to ~ 65% for the most luminous ones (Mt x ^ —24, 
bin 1). This gradual reduction of the fraction of stellar mass 
contributed by the stars for the more luminous galaxies oc- 
curs at expenses of the accretion of the stellar component of 
merging satellites ('sat' contribution), indicating that merg- 
ers become important for the more massive galaxies. 

We now analyse the evolution with redshift of the ac- 
cumulated stellar mass contributions from the different pro- 
cesses, within different magnitudes bins. The evolution of 
the corresponding mean fractions of stellar mass are shown 
in Fig. [4] for galaxies with magnitudes corresponding to bins 
1, 3 and 6; these fractions are given with respect to the total 
contribution to the stellar mass of galaxies at the redshifts 
considered. Although we have taken into account both the 
'stars' and 'sat' components arising from major dry mergers, 
we found that the contribution of the former to the whole 
stellar mass of galaxies is negligible. Hence, in the following, 
we simply refer to this process as 'major dry merger', mean- 
ing that only the stellar mass from the 'sat' component is 
being considered. 

From a first inspection of these plots, we see that qui- 
escent SF appears as the dominant process that contributes 
to the stellar mass of galaxies within all magnitude bins and 
at all redshifts. However, we can see that the corresponding 
stellar mass fractions decrease monotonically with redshift 
as the cold gas reservoir in each galaxy is exhausted. 



The relative contribution of the rest of the processes, 
broadly grouped in disc instabilities and mergers, change 
considerably as we move along the RS. The left panel of 
Fig. [4] shows the mean stellar mass fractions for galaxies 
with low luminosity (—18 < Mt x ^ —17, bin 6). We can 
clearly see that the second larger contribution to the stellar 
mass of these galaxies comes from the stars formed during 
minor wet mergers (i.e. the 'stars' component); this mass 
fraction increases with redshift reaching a value of w 25 per 
cent at z fa that has already been achieved at z ~ 2, 
indicating that minor wet mergers do not play a significant 
role since that epoch. These fractions are followed by those 
corresponding to the contribution of the major wet merger. 
As we will see, the contribution of this particular process is 
very small or almost negligible for more luminous galaxies, 
therefore, for the sake of clarity, we do not split this process 
in the two different components. The mass contributions 
provided by disc instabilities and other kind of mergers are 
negligible. 

As we move up along the RS to higher luminosities, the 
situation changes appreciably. The mass fractions for galax- 
ies in the critical magnitude bin 3, where the detachment of 
the bright end of the RS from the general linear fit occurs, 
are shown in the middle panel of Fig. [3] At these magni- 
tudes, the relevant contribution to the stellar mass content 
(apart from quiescent SF) arises from minor dry mergers, 
while the relevance of the 'stars' component in minor wet 
mergers decreases considerably from z fa 4, becoming negli- 
gible at z = 0, in contrast to the behaviuor manifested by 
this process for galaxies in bin 6. We also find a mild contri- 
bution of major wet mergers at low redshift. All these merger 
events are responsible of ~ 20% of the stellar mass received 
by galaxies with magnitudes in the range —21 < Mt x ^ —20 
at z — 0. 

For the most luminous galaxies (—24 < Mt± ^ —22, bin 
1), the trend of the mass contribution of minor dry mergers 
is reinforced with respect to the one shown for bin 3, gain- 
ing importance from earlier epochs (i « 6), and reaching a 
higher value of ~ 22% at z = 0. The contribution of ma- 
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Figure 5. Number of mergers suffered by the galaxies in a C15 
cluster through their whole lifetime as a function of the baryonic 
mass ratio of mergers; different curves correspond to magnitude 
bins 1, 3 and 6, which contain different number of galaxies. 



jor dry mergers also increase notably since z sa 2 giving a 
fraction of « 10% at Z — 0. Thus, the whole contribution 
of these two relevant processes at the present epoch leads 
to a fraction of sa 32%, consistent with the result shown in 
Fig. E] 

Although the previous description is based on the mean 
mass fractions of galaxies in the massive C15 clusters, which 
clearly show the change in slope of the RS, similar results are 
obtained for the small simulated clusters C14 (not shown). 
In the latter, however, statistics become very poor at high 
luminosities, with only a few objects, if any, in the highest 
magnitude bin. This g eneral behaviour is in agreement with 
the scenario infered bv lMarti'nez et al.l (|2010r ) from the anal- 
ysis of galaxies in groups and in the core of X-ray clusters, in 
which these galaxies populate nearly the same CMR. Their 
result supports previuos works that indicate th at the CMR 
is independent of the environm ental overdensity IjHogg et al.l 
12004 lLopez-Cruz et alj|2004h . 



4.2 Mass ratio of mergers resolved in the model 

The contribution of minor mergers to the stellar mass is 
relevant for galaxies in all magnitude bins, as it is evident 
from Fig. [4] The mass ratio of minor mergers that are re- 
solved by the model depends on the mass of the galaxies 
involved. Given the mass resolution for DM haloes in the N- 
body/SPH simulations used, a very conservative upper limit 
for the mass resolution in galaxies would be ~10 9 /i _1 M© 
(the minimum resolved halo mass times the baryonic frac- 
tion), but due to the complex interplay of processes such 
as gas cooling, feedback and strangulation, we find that 
our model produces objects with stellar mass as small as 
1.5 x 10 s h~ x Mq. Therefore, the minimum mass ratio re- 
solved in the model (corresponding to the merger of such an 



object with a BCG, with stellar mass of ~ 1.5 x 10 12 /i _1 M 
in our simulations) is ~ 10 -7 . However, it is likely that our 
model is incomplete for such small objects. 

Comparing the results of the model with e.g. 
the observed lumino sity function for galaxy clusters 
IjDePropris et al.l 120031 ). it seems our model is complete 
for galaxies with absolute magnitude < —15 in the bj- 
band. This corresponds to galaxies with a mean stellar mass 
~ 2 x Kfh" 1 Mq. Consequently, for the case of the BCGs, 
we expect our model to be complete in mergers with mass 
ratio >10 -4 . For galaxies in the C15 clusters, the mini- 
mum mass ratios resolved in the model are of the order 
of 4 x 10" 5 , 10" 4 and 6 x 10" 3 for bins 1, 3 and 6, re- 
spectively. The ratios indeed increase as one considers less 
luminous galaxies. When selecting the most luminous galax- 
ies in bin 1 in the C14 clusters, the minimum mass ratio is 
higher but of the same order of magnitude, ~ 2.5 x 10 -5 . 

In general, the number of mergers suffered by galaxies 
along their whole lifetime increases rather monotonically as 
mergers of lower mass ratio are considered, as can be seen 
from Fig. [S] where the results are shown for magnitude bins 
1, 3 and 6. Minor mergers (defined as those with mass ratio 
lower than 0.3; Section I2.3|l are dominant for all galaxies 
along the RS. For the least luminous bin, the maximum is 
not reached for the least value of the mass ratio, as in the 
more luminous ones, reflecting the influence of the limit in 
the mass resolution of the simulations. Each of these bins 
contains different number of galaxies, as indicated in the 
plot, being larger for the least luminous one. Taking into 
account this aspect, the average number of mergers is larger 
for brighter galaxies. 



4.3 Influence of gas fraction thresholds 

From the analysis of the stellar mass fractions discussed in 
Subsection 14.11 it is clear that the evolution of galaxies in 
the bright end of the RS (those brighter than the critical 
value (Mt-l ~ —20) is driven by the dry contribution of 
both minor and major mergers. On the other hand, fainter 
galaxies are affected by the wet component of mergers. 

These conclusions are in agreement with those obtained 
by ISkelton et al.l l|2009r ) , who find that most of the galaxies 
in the bright end of the CMR of local early-type galaxies 
have undergone dry mergers, in contrast to faint galaxies. 
Their results are based on a simplified model that isolates 
the effects of the dry merging in the colours of the RS galax- 
ies. The merger histories and gas fraction s of galaxies are 
extrac ted from the semi-analytic model of lSomerville et al.l 
(2008). Their model is able to reproduce the change of 
slope in the bright end of th e CMR as seen in local early- 
type galaxies from the SDSS. ISkelton et al.l affirm that this 
change in slope, and the magnitude at which the break oc- 
curs, depend strongly on the assumption of the gas frac- 
tion threshold below which mergers are considered dry. This 
threshold is allowed to vary between 10 and 30 per cent. 

In our model, on the contrary, the tilt of the RS to- 
wards bluer colours has only a weak dependence on the 
gas thresholds adopted to distinguish between wet and dry 
mergers (defined in Section EOI as / gas , minor and / gas , major 
for the minor and major mergers, respectively). The pa- 
rameter /g as , major has no influence on the effect of major 
mergers, because such mergers always lead to the forma- 



12 Noelia Jimenez et al. 



Magnitude bins 
2 3 4 5 



1.0 



0.5 



0.0 



-0.5 



-1.0 



-24 < M T1 < -22 

-22 < M T1 < -21 

-21 <M T1 <-20_ 

... -20<M T1 <-19 

_ -19<M T1 <-18 

„, -18<M T1 <-17 

... -17<M T , <-16 




0.30 



0.25 



E 0.20 



0.15 



0.10- 



0.05- 



0.00 



-24 



- 








I I I . 


- 


since z= 






accreted stellar mass 


" 








('sat' contribution) 


- 








stellar mass formed _ 


- 


;, : 






('stars' contribution) - 


- 




''• 


'• 


• 


— 








+ * — 


- 


..---' 

» 


.•* 


* - 


* * 


- 








m. 


- 








i I i I i '"1 * 



-23 



-22 -21 -20 -19 -16 
Absolute magnitude 



-17 



Figure 6. Evolution with redshift of iron abundance of the stellar 
component of galaxies that at z = lie within a given range of 
magnitudes; metallicity values are a result of the complex com- 
bination of all processes included in our model SAG. Averages 
at each redshift have been estimated considering galaxies in C15 
clusters. Different lines represent the mean values for galaxies 
within different magnitude bins. 



tion of a spheroid and the consumption of the cold gas in 
a starburst, regardless of the amount of gas available in the 
merging galaxies. The threshold / gas , major only has the pur- 
pose of classifying major mergers as either dry or wet, and 
thus only affects the percentages of the 'stars' and 'sat' com- 
ponents contributed by these processes to the stellar mass 
of galaxies. We have verified that changing this threshold 
within the range 0.2 ^ / gas , major ^ 0.6 results in no signifi- 
cant changes with respect to the results presented in Figure 



[4] obtained for / g 



= 0.4. On the other hand, the 



value /gas, minor = 0.6 for the minor mergers has been tuned 
together with the rest of the free parameters of the model 
in order to reproduce several observed galaxy properties, as 
described in Section [2, 3 1 

We would like to emphasize that our model was not con- 
structed to produce a break in the RS; this arises naturally 
as a consequence of the physical processes included in the 



model. In particular, varying the value / ga 



over the 



aforementioned range, a break in the bright end of our sim- 
ulated RS is always present, without any appreciable change 
in the values of the slope or of the critical bre ak magnitude. 
In this respect, our findings differ from those of lSkelton et al.l 
(2009), who claim to see a strong dependence of the break 
on the adopted gas threshold; moreover, we do not see a sig- 
nificant dependence from their results (see their Fig . 1). In 
any case, both our results and those of lSkelton et al.l (2009) 
strongly support the fact that the dry contribution of the 
merger processes are relevant in the formation of the bright- 
est galaxies in the RS, both in field and cluster environ- 
ments. It has been suggested that dry mergers are expected 
to only increase the mass of the remnant without changing 
their colours because there is no associated SF, thus caus- 
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Figure 7. Top panel: stellar mass fractions of 'stars' and 'sat' 
contributions accumulated since z = 1 for each magnitude bin, 
normalized with the total stellar mass of galaxies at z = within 
a given bin; galaxies in the C15 clusters are considered. Bottom 
panel: Average iron abundance of galaxies at z = as given by the 
complete evolution traced by the semi-analytic model (solid line), 
for galaxies within different magnitude bins, compared with the 
metallicity of the mass contributed by the stars formed since z = 1 
(dashed line) and of the stars already present in the satellites 
accreted since z = 1 (dot-dashed line). 



ing th e detachment of the bright end from the general linear 
trend (JBernardi et aLhOOTUSkelton et aJBoQJ ). In what fol- 
lows, we explain this claim by means of our detailed model 
of galaxy formation, focusing on cluster galaxies. 
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4.4 Metallicity evolution 

The evolution of the mean value of the stellar iron abun- 
dance of galaxies that at z — lie within a given range 
of magnitudes is shown in Fig. [5] These metallicities are 
the result of the complete combination of all the processes 
considered in the semi-analytic code; averages have been es- 
timated for each magnitude bin considering the three C15 
clusters. We can see that, regardless of the galaxy magni- 
tude, the mean metallicity of galaxies increases smoothly 
with redshift, being always higher for the more luminous 
galaxies. In particular, at z — 0, galaxies in the bright end 
(Mt x > —20) reach metallicity values within a narrower 
range (tzs 0.15 dex) than the rest of the galaxy population. 
This fact is directly linked with the similar colours that char- 
acterize these galaxies making them depart form the general 
trend of the RS. These chemical abundances seem to be in 
place since z = 1. This is supported by observations of high 
redshift clusters, which show that the slope and scatter in 
the CMR for morphologically selected ET Gs show little or 
no evidence of evolution out to z » 1.2 fBla keslee et al.l 
l2003l : lMei et alJl200d : IJaff et alJl2010h . 



4.5 Mass assembly and merger rates since z — 1 

Taking into account that metallicity values of galaxies in the 
bright end have already been achieved at z — 1 (Fig. [6]), we 
consider the accumulated mass contributed from both com- 
ponents ('stars' and 'sat') in the redshift range < z < 1. 
The relative contribution of these two components are shown 
in the top panel of Fig. [7] for galaxies in different magnitude 
bins. These mass fractions are normalized with the total 
stellar mass of galaxies at z = within a given magnitude 
bin, because we are interested on the relative importance of 
the mass accumulated since z — 1 to the stellar mass of the 
galaxy at the present epoch. Note that this stellar mass is 
a result of the whole evolution traced by the semi-analytic 
code, including its reduction as a result of the mass recycling 
arising from stellar evolution. However, by construction, the 
variables used to analyse the mass fractions shown in Fig. 
[7] (as well as in Figs. [3] and [4]) save the information of the 
mass of stars borned in a galaxy or added to it during merger 
events, but they are not updated to take into account the 
recycled mass. This fact specially affects the stellar mass 
fractions contributed by the mass of new born stars ('stars' 
component) since z — 1, which are shown by the dashed 
line in the top panel of Fig. [Jj As we can see, these fraction 
ranges from ~ 0.05 for the brightest galaxies to w 0.15 for 
the faintest ones. However, they are in fact even lower taking 
into account that they are normalized with the total stellar 
mass of galaxies at Z — 0, which have lost « 35 per cent 
of their mass as a result of the evolution of stars belonging 
to the high mass tail of a Salpeter IMF (which is the one 
adopted in our semi-analytic model). 

These fractions of the stellar mass formed in situ since 
z — 1, in which all physical procesesses contributing with 
cold gas from which new born stars originate are taken 
into account, can b e compared with estimations done by 
iKavirai et al.l (|2008r ). They compute the recent star forma- 
tion (RSF; mass fraction of stars that form in the galaxy 
within the last 1 Gyr of look-back time in its rest frame) in 
the high redshift (0.5 < z < 0.1) early-type population from 



rest-frame UV and optical colours and com bine their re- 
sults with similar studies in the local universe IjKavirai et al.l 
120071 ). thus being able to estimate the star formation history 
of the early- type population between z — and 1. For the 
range 0.5 < z < 1, early type galaxies belonging to the RS 
(NUV — r > 4) typically show RSF values less than 5 per 
cent and the reddest early-types (which are also the most 
luminous) are virtually quiescent with RSF values of ~ 1 per 
cent. 

On the other hand, IKavirai et al.l J2007h find that the 
SF activity in the most active ETGs (those in the blue peak 
of the colour-magnitude distributions) has halved from an 
average RSF of ~ 11 per cent at z = 0.7 to ~ 6 per cent at 
the present day. From these results, they infer that 10 — 15 
per cent of the stellar mass in luminous early-type galaxies 
(—23 < Mv < —20.5) may form after z = 1, while less lumi- 
nous galaxies may form 30 — 60 per cent of their mass within 
this time-scale. The fractions of the stellar mass formed in 
situ since z — 1 given by our model (dashed line in Fig. [7J 
are a factor of 2 smaller than these values, within ranges of 
~ 0.05 — 0.09 and ~ 0.09 — 0.15 per cent for galaxies in the 
bright and faint end of the RS, respectively. This difference is 
expected, taking into account that we are analysing cluster 
galaxies where the peak of the star formati on takes place 
at ear lier times, while those in the study of IKavirai et al.l 
(2008) belong to a variety of environments. Thus, our re- 
sults are in very good agreement with these observational 
estimates. 

The stellar mass fractions contributed by the accreted 
stellar mass during mergers are not so affected by the defi- 
nition of the accumulated quantities, since the stellar popu- 
lation of satellites have already suffered from the recycling 
processes by the time they are accreted, assuming that the 
bulk of their stars where formed at z > 1. Indeed, if a satel- 
lite galaxy accreted within the redshift range < z mcrgcI < 1 
is created as early as z — 7, then, by the time it merges 
with the central galaxy it has already lost ~ 36 per cent of 
its original mass retaining only the low mass (slowly evolv- 
ing) stars. Then, the stars that will be added to the central 
galaxy will evolve very slowly and more than 93% of those 
stars will still be alive at z = 0. Even in the case of a rela- 
tively young galaxy created at z — 2 and accreted during the 
time interval < z mcrgor < 1, we find that more than 90% 
of the stellar mass added during the merger is still alive at 
z = 0. Keeping these caveats in mind, we can conclude from 
Fig. [7] (upper panel) that ~ 10 — 20 per cent of the mass 
of galaxies in the bright end of the RS arises from satellites 
accreted since z — 1, with very few stars being formed in 
situ. 

Observational quantification of the impact of mergers 
in the mass assembly of early type galaxies yields to larger 
increments of the mass of galaxies since z = 1 with respect 
to what we have found. However, we have to keep in mind 
that our results correspond to cluster galaxies, characterized 
by a faster evolution , as al ready mentioned. In particular, 
lLopez-Saniuan et at.l (J2010j) find that the brightest early- 
type galaxies increase their mass by about 42 ± 8 per cent as 
a result of mergers since z — 1. They use kinematically con- 
firmed close pairs from the VIMOS VLT deep spectroscopic 
redshift survey to study the role of minor mergers to mass 
assembly of normal Lb > L* B galaxies, which can belong 
to different environments. The fact that the mass fraction 
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received by massive early-type galaxies in their sample is a 
factor of two larger than the one estimated for our brightest 
galaxies is a consequence of the lower number of mergers 
that cluster galaxies suffer since 2 = 1. 

The estimation of the ave rage number of mergers per 
red galaxies since z — 1 done bv lLopez-Saniuan et at.l (|2010l ) 
is 1.3 ± 0.3, that is splitted in 0.7 ± 0.1 major mergers and 
0.6 ± 0.2 minor ones. Their criterion to separate between 
minor and major mergers is based on the ratio of B-band 
luminosity of the merging components which is used as a 
proxy of mass; minor mergers are characterized by the lu- 
minosity ratio ranging from 0.1 to 0.25. In our model, minor 
mergers are defined as those having mass ratios less than 
0.3, whithout impossing a lower limit. Thus, we find that 
galaxies in the bright end of the RS (Mt± ^ —20) have suf- 
fered an average number of mergers of ~ 4.88 since 2 = 1, 
which is almost four t i mes l arger that the value obtained by 
lLopez-Saniuan et at.l (|201(J ) . Therefore, we repeat our esti- 
mation restricting to minor mergers with mass ratios larger 
that 0.1. Thus, we find that the average number of mergers 
suffered by galaxies in the bright end of the RS is consid- 
erably reduced, with ~ 0.33 minor mergers per galaxy. The 
value is a factor of t wo sm aller that the one inferred by 
lLopez-Saniuan et at.l (|2010r i. but again, we have to keep in 
mind that our sample contains cluster galaxies. 

The cumulative fraction of stellar mass contributed by 
minor mergers since 2 = 1 that we compute, ranges from 
~ 0.12 to 0.07, for the most and the least luminous galax- 
ies in the bright end of the RS, respectively, which are in 
very good agreement with the estimation of 0.1 done by 
lLopez-Saniuan et at.l (|2010P l. These values have been ob- 
tained in the same way as those shown in the upper panel 
of Fig. [7] but restricting to the contribution of the 'sat' 
component of minor mergers, including those with mass 
ratio less than 0.1. This indicates that the mass contri- 
bution of minor mergers with such a low mass ratio, al- 
tough being highly numerous, contribute with a very small 
percentage to the stellar mass of the central galaxy. On 
the other hand, the average number of major merger suf- 
fered by the brightest galaxies in the RS since z — 1 is 
~ 0.19, a factor of three sm aller than the one obtained by 
lLopez-Saniuan et at.l (|2010| ). which explains that the mass 
fraction contributed by all kind of mergers to the brightest 
galaxies in the RS (~ 0.2) is a factor of two smaller than 
observed bv lLopez-Saniuan et at.l (|2010j) (~ 0.42). However, 
this value of the number of major mergers is in better agree- 
ment with other empirically determined major merger rates, 
which indicates that major merger activity in galaxies with 
masses comparable with the ones in this paper, becomes in- 
frequent after 2 = 1 (e. g. IConselice et aUboOa ; iBundv et al.1 
l2009l ; |jogee et alj|2009h . 

We also compare our results with those obtained by 
iNierenberg et al.1 (|2011i ) . They study the spatial distribution 
of faint satellites of early-type galaxies at intermediate red- 
shifts (0.1 < z < 0.8) selected from the GOODS field. Satel- 
lites are detected from high resolution HST images consider- 
ing those that are up to 5.5 magnitudes fainter than the host 
galaxy and as close as 0"5/2.5 kpc to the host. The major- 
ity of their host galaxies are massive elliptical galaxies with 
typical stellar masses of ~ 1O 1O ' 5 M0 that have on average 
l-Tto'ij satellites. This inferred satellite number can easily 
account for the number of minor mergers (0.4 — 0.5) that 



60 per cent of their host galaxies are expected to undergo 
in the time span of the study. This number of mergers has 
been obtained by using the fitting form for the mean merger 
rate per halo as a function of halo mass, progenit or mass ra- 
tio and redshift derived bv lFakhouri et al.1 (|2010P ) from halo 
merger statistics. Although a bit lower than these infered 
values, the result from our model (~ 0.33 minor mergers 
per galaxies since 2 = 1) is consistent with them. 



4.6 Metal contribution since 2 = 1 

An interesting aspect of the plot in the lower panel of Fig. [7] 
is that the contributions of the 'stars' and 'sat' components 
cross each other at the break magnitude Mt x = —20, clearly 
indicating that the stellar population accreted during late 
mergers determines the properties of galaxies in the bright 
end of the RS. From the discussion related to this plot, we 
can safely assume that more than 90% of the stellar mass 
formed and added by mergers since 2 = 1 is still alive at 
2 = 0, determining the metallicity and colour of galaxies in 
the RS at the present epoch. We estimate the metallicity 
of the stellar mass contributed since z = 1 by the these 
two components, [Fe/H] star s and [Fe/H] sa t, considering the 
mass of iron, Fe^° cc ° mp , and hydrogen, H^° cc ° mp , defined 
in Section[4] The average value of these metallicities for C15 
cluster galaxies within different magnitude bins are shown 
in the lower panel of Fig. [7] They are compared with the 
average metallicity achieved by the stellar mass of galaxies 
at z = according to the whole evolution traced by the semi- 
analytic code (solid line), which lies in the range —0.65 < 
[Fe/H] total % 0.05 for the least to the most luminous galaxies 
in the RS. The metallicity of the new stars formed since z = 
1 (dashed line) ranges from [Fe/H] sta r ~ —0.45 for the least 
luminous galaxies to ~ 0.5. In particular, for galaxies in the 
bright end of the RS, with magnitudes —24 < M Tl < -20 
(bins 1, 2 and 3), the metallicity of the 'stars' component 
is above solar (0.25 < [Fe/H] st ars < 0.5) as the result of 
the highly chemically enriched cold gas available for star 
formation in more massive galaxies. However, as we have 
seen, the mass fraction contributed by processes involving 
star formation is very low for these galaxies (see top panel 
of Fig. [7} , therefore, these metal contributions hardly affect 
the metallicity of the whole galaxy. On the other hand, the 
'sat' component of satellites accreted since 2 = 1 represents 
an appreciable fraction of the galaxy mass at z = for 
galaxies in the bright end (sa 10 — 20 per cent), but the values 
[Fe/H] sa t (dashed-dotted line in bottom panel of Fig. [7} are 
subsolar for galaxies in all magnitude bins, being w 0.15 
dex lower than the corresponding galaxy metallicity. Hence, 
dry mergers help increase the mass of galaxies in the bright 
end of the RS without considerably affecting their average 
metallicities. This upper limit in the metallicity achieved by 
galaxies in the bright end fixes their colour, being bluer than 
expected if star formation from highly chemically enriched 
gas were relevant in their final evolution, thus explaining the 
departure of the bright end of the RS from a linear fit. 



4.7 Discussion 

Our findings related to the low metallicity of the stellar mass 
of satellites accreted al low redshifts (z < 1), with respect to 



The bright end of the CMR 15 



the chemical abundances found in the r emnan t galaxy, are 
in line with the results of lLagos et al.l (120091 ) . They show 
a clear correlation between the metallicity of galaxies and 
the mass of their host DM haloes, with central galaxies be- 
ing more chemically enriched than their surviving satellites. 
This trend is a direct consequence of the hierarchical clus- 
tering process, where more massive galaxies sit in initially 
higher density peaks and, as a consequence, start to form 
earlier on average, allowing them to acquire higher amounts 
of metals. 

Our results also support the scenario proposed by 
iMartmez et al.l IJ2010I ) for the formation of the CMR in mas- 
sive systems. Their analysis of observational results indicate 
that the CMR is populated by red ETC that formed the 
bulk of their stars during the early stages of massive halo 
assembly, on one hand, and red galaxies that passed most of 
their lives inhabiting poor groups or the field and fell into 
massive systems at lower redshifts, on the other. Our study 
suggests that part of the later population are the satellites 
intervening in the dry mergers processes, that play a signif- 
icant role in the late evolution of massive galaxies. 

The fact that massive cluster galaxies are fed by mi- 
nor dry mergers whereas fainter galaxies along the CMR 
are mainly affected by minor wet ones is due to the in- 
fluence of the environment in which the host galaxy and 
their satellites reside. More massive galaxies tend to inhabit 
denser regions, where it is highly likely that infalling satel- 
lites are stripped of their cold gas before the merger as a 
result of environmental effects. Hence, the modelling of the 
environmental processes that affect the cold gas content of 
galaxies could influence our conclusions regarding the frac- 
tion of mass received by galaxies from dry or wet mergers. 
In our model, as in most previous semi-analytic modelling, 
the hot gas haloes of galaxies are completely and immedi- 
ately stripped as soon as a galaxy becomes a satellite of a 
group or cluster. The affected galaxy cannot replenish its 
cold gas any further, and with a typical SFR it exhausts 
its co ld gas in a few Gyr (|Gallagher. Hunter fc Bushouse 
1989). This process is called 'strangulation' (jLarson et al. 
19801 ). and provides an e xplanation for the SFR gradients 



in clusters lIBalogh et al.ll2000l) and the c olour evolution of 
cluster ga laxies IlKodama fc BowerlbOOli ). However, recent 
work (e.g. IWeinmann et al.ll2006r ) has shown that, compared 
to SDSS observations, semi-analytic models produce higher 
fractions of red satellites. This is likely due to the crude mod- 
elling of the removal of the hot gas haloes. If the hot gas is 
removed more gradually, satellite galaxies can prolong their 
star formation resulting in bluer colours. It is not known 
which physical process is actually responsible for the removal 
of the hot gas, but one lik ely candidate is ram-pre ssur e strip- 
ping ( RPS), as shown bv lMcCarthv et al.1 |200Sl) and lBekkil 
(I2009F ). 

This satellite overquenching problem raises the issue 
that perhaps the semi-analytic model used overpredicts the 
number of dry mergers, as satellite galax ies in the model 
would be too gas-poor. I Font et al.l i|2008l ) has shown that 
considering a gradual removal of hot gas by RPS, a better 
match to the observed fractions of red galaxies is obtained. 
One would conclude that with such a modification, the num- 
ber of wet mergers should increase; however, the situation is 
not so simple. If the gas content of satellite galaxies is reg- 
ulated by RPS, the colour of the satellite will be bluer as a 



result of a more prolonged star formation, but this does not 
necessarily mean that a higher mass of gas is contributed 
to the mergers. In the central regions of galaxy clusters the 
ram-pressure exerted by the ICM is stron g enough to signif - 
icantly affect the gas content of galaxies (|Tecce et al.ll2010r ) 
and by the time the satellite merges with the central, not 
only its hot gas but also the cold gas could be stripped 
away by RPS, becoming part of the ICM. Thus, the gas 
content of the satellite would be mostly removed prior to 
the merger (|Sofueill994l ; lGrebel. Gallagher fc Harbeckll2003l ) 
and the number of wet mergers might not inc rease signifi- 
cantly . We have updated the model described in lTecce et al.l 
(2010) to consider the RPS of the hot gas, and we intend to 
explore this issue further in a forthcoming paper. 

Finally, it is interesting to note that the break magni- 
tude in which there is a change of the slope in the RS {Mt 1 ~ 
—20), corresponds to a galaxy mass of M* ~ 10 10 ft -1 Mq, 
also marks the difference between two regimes in the depen- 
de nce of metallic i ty gra dients with galaxy mass, as found 
bv lSpolaor et al.l (|20091 ). Their results suggest that cluster 
galaxies above the mentioned mass threshold might have 
formed initially by mergers of gas rich disc galaxies and then 
subsequently evolved via dry merger events, in concordance 
with our findings. 



5 SUMMARY AND CONCLUSIONS 

In this work we have studied the physical processes respon- 
sible for the detachment of the bright end of the RS of clus- 
ter galaxies from the linear fit determined by less luminous 
galaxies. The RS is constructed by selecting galaxies accord- 
ing to the empir ical redshift-dependent criterion given by 
iBell et al.l (J200J) to separate the blue and red sequences 
in the colour-magnitude diagram. We use the s emi- analytic 
model of galaxy format ion and evolution SAG (jLagos et al.l 
120081 ; iTecce et al.l I2010T ) . combined with cosmological non- 
radiative A^-Body/SPH simulations of galaxy clusters in a 
concordance A Cold Dark Matter universe. We consider 
two sets of simulated clusters with virial masses in the 
ranges ~ (1.1 - 1.2) x 10 14 /i~ 1 M Q (C14 clusters) and 
~ (1.3 - 2.3) x 10 15 ft" 1 M Q (C15 clusters). 

We find a very good agreement between the general 
trends of the simulated RS and observed CMR of ETGs in 
four different colour-magnitude planes. The simulated rela- 
tion presents a break at approximately the same magnitude 
in the different systems (M^ roak ~ M^ roak ~ M^ roak « -20, 
galaxy mass of ~ 10 10 Mq), with more luminous galaxies 
having almost constant colours (Washington C — T\ ~ 1.5) 
(see Fig. [TJ instead of following the linear trend infered from 
the less luminous ones. 

This aspect is consistent with a non-linear fit found 
from different set of observation s (Baldrv et al. 2006; 
iJanz fc Liskerl 120091 ; ISkelton et all 120091 ). The change of 
slope in the bright end of the RS emerges naturally from our 
model once it has been calibrated to satisfy several obser- 
vational constraints simultaneously. One of the obervational 
relationships that is closely followed by the simulated galax- 
ies is the luminosity- metallicity relation (see Fig. [SJ. This 
fact supports the use of the chemical history of galaxies as a 
tool to help understand the development of the RS and its 
break at the bright end. 
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Thus, we focus on the analysis of the mass and metals 
contributed to each galaxy by different physical processes: 
quiescent star formation, and starbursts during disc insta- 
bility events and mergers. We study the mean properties 
of galaxies in seven different magnitude bins; galaxies are 
grouped according to their Mt x magnitudes. From the anal- 
ysis of the accumulated stellar mass contributed by different 
processes, we obtain information about the way in which 
their relative importance change with redshift, for galaxies 
in a given magnitude bin, and the dependence of their rele- 
vance with galaxy luminosity. We find that: 

• Quiescent SF appears as the dominant process that con- 
tributes to the stellar mass of galaxies within all magnitude 
bins and at all redshifts. The corresponding stellar mass frac- 
tions decrease monotonically with redshift as the cold gas 
reservoir in each galaxy is exhausted. 

• Galaxies with low luminosity ( — 18 < Mt x ^ —17, bin 
6) receive a large contribution from the stars formed during 
minor and major wet mergers, with the former processes 
giving place to a larger mass fraction, which increases with 
redshift and reaches a value of ~ 25 per cent at z « 0. This 
value has already been achieved at z fa 2, indicating that 
minor wet mergers do not play a significant role since that 
epoch. 

• The most luminous galaxies (—24 < Mt x Sj —22, bin 1) 
have a considerable contribution from minor dry mergers, 
which become important since z ~ 6; the corresponding 
mass fraction reaches a value of ~ 22 per cent at z = 0. 
The contribution of major dry mergers is also relevant for 
these galaxies since z » 2, giving a fraction of « 10 per cent 
at z = 0. 

• The transition between the relative importance of the 
contribution of wet and dry mergers occurs for galaxies 
within the range of magnitudes (—21 < Mt x ^ —20, bin 3; 
M^ ~ 1O 1O M0), that is, galaxies characterized by the break 
magnitude, where the detachment towards bluer colours in 
the RS takes place. At these magnitudes, minor dry merg- 
ers start to gain importance, while the relevance of the stars 
formed during starbursts in minor wet mergers decreases 
considerably from z ~ 4, becoming negligible at z — 0, in 
contrast to the behaviuor manifested by this process for the 
least luminous galaxies (bin 6). We also find a mild con- 
tribution of major wet mergers at low redshifts. All these 
merger events are responsible of « 20 per cent of the total 
mass received by galaxies. 

From the analysis of the relative importance of the stel- 
lar mass fractions contributed by different processes, it is 
clear that the evolution of galaxies at the bright end of the 
RS is driven by the dry contribution of both minor and ma- 
jor mergers, while fainter galaxies are affected by the wet 
component of mergers. This situation is already established 
at z ~ 1. The metallicities of galaxies in the bright end do 
not change since that redshift, although a mild increase is 
detected for lower luminous galaxies where star formation 
from chemically enriched gas keeps taking place. 

Considering the mass accumulated since z = 1 from 
both the stars formed in situ and and the stellar mass al- 
ready present in the accreted satelites, we find that the latter 
represents a fraction of the present stellar mass of galaxies 
in the bright end of the RS that ranges between w 10 and 
20 per cent. This fraction is larger than the fraction of stars 



formed during recent star formation. Altough the metallicity 
of the stars formed since 2«lis high (« 0.2 — 0.5 dex above 
the mean iron abundances of galaxies at z = 0, see Fig. [7] 
bottom panel), their mass contribution is so low (~ 0.05) 
that they do not affect the mean metallicity of galaxies. On 
the contrary, dry mergers provide a higher fraction of stellar 
mass than that generated by recent SF but with low content 
of metals (~ 0.2 dex below the mean metallicity of galax- 
ies) thus, their contribution only increase the stellar mass of 
galaxies without changing their mean metallicities. Hence, 
galaxies in the bright end, which have negligible star forma- 
tion activity since z f» 1, reach an upper limit in their abun- 
dances that in turn fixes their colours, taking into account 
the small spread in age of these systems. The departure of 
the bright end of the RS from a linear fit can be explained 
by the effect of the contribution of minor and major dry 
mergers, which increase the mass/luminosity of the galaxies 
but without changing their colours. 

These conclusions are valid for both high and low mass 
galaxy clusters. We have demonstrated that the local prop- 
erties of galaxies on the RS seem to be established since 
2 ~ 1. A detailed analysis of the features of the RS, such 
as its dependence with environment and redshift, will be 
addressed in a series of forthcoming papers. 
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